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The model is used to describe low concentration of bivalent ions adsorp- 
tion in the systems of environmental importance. A theory has been de- 
veloped based on a picture of a heterogeneous solid surface, with different 
adsorption sites — the outermost surface oxygens. The theoretical study in- 
dicates that the binding-to-surface energies of different surface complexes, 
formed on different surface oxygens, vary in an independent way from one 
surface oxygen to another. Furthermore, the idea of the Four Layer Model, 
where a new layer was reserved for the bivalent metal ions, was accepted. 
Cations and anions of basic electrolyte were placed in the same layer as in 
the popular Triple Layer Model. 


1. INTRODUCTION 


The models of surface complexation are the most frequently applied ones to 
describe ion adsorption within the double electrical layer formed at metal ox- 
ide/electrolyte interfaces. 

The adsorption of Me** metal ions at electrolyte/oxide interfaces at low ion 
concentrations is a subject of rapidly growing interest. These are problems of 
great practical importance: the adsorption of highly poisoning cations of some 
heavy elements in soil; the adsorption of radioactive ions in soil in the areas 
where nuclear plants are located; the adsorption of radioactive ions on corroded 
parts of nuclear plant installations. In all these systems, concentrations of the 
ions are low, and their adsorption characteristics are strongly affected by the 
oxide surface energetic heterogeneity. 

Of all the models published so far, the most frequently used was the so 
called the "Triple Layer Model" (TLM) [1], which according to the up-to-date, 
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more realistic outlook has been modified by us by assuming that the surfaces of 
real oxides are geometrically deformed and energetically heterogeneous [2-5]. 

Our previous paper [5] presented the theoretical studies based on the TLM 
model in which the bivalent metal ions formed so called 'inner-sphere' com- 
plexes in the surface layer, where the potential determining ions H* are ad- 
sorbed. 

Such an assumption was made in the studies by Hayes and Leckie [6-9]. 
They stated that it is more consistent with reality than assuming formation of so 
called 'outer-sphere' complexes in which the bivalent metal ions are adsorbed in 
the same layer as cations and anions of the basic electrolyte. Although the 
model presented in paper [5] was more realistic, it allocates ions of very differ- 
ent sizes and degrees of hydration in the same plane (especially protons and 
bivalent ions). However, introducing another plane would involve a large in- 
crease in mathematical complexity. 
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Figure 1. Diagrammatic presentation of the model of the metal oxide/electrolyte inter- 
face. Wo, 5p - the surface potential and the surface charge density in the O-plane; Wm, Ôn - 
the potential and the charge coming from the bivalent cations adsorbed; ye, 6p - the 
potential and the charge coming from the specifically adsorbed ions (cations C* and 
anions A`) of the basic electrolyte; wy, 5, - the diffuse layer potential and its charge; c1, 
c2- the electrical capacitances, constant in the regions between planes 


A four layer model for bivalent cations adsorption... 189 


So, recently we have accepted the idea of the Four Layer Model (FLM) in- 
troduced to literature by Bowden et al. [10,11] as well as in the papers by Bar- 
row [12,13] and Sposito [14]. A new layer, (the fourth one as the name indi- 
cates, but situated as the second, next to the surface layer "0" where protons are 
adsorbed), was reserved first for the bivalent metal ions. Cations and anions of 
basic electrolyte were placed in the same layer as in the Triple Layer Model. A 
schematic picture of FLM taking into account the bivalent cation adsorption is 
shown in Figure 1. That model is used here by us to describe low-concentration 
of bivalent ions in systems of environmental importance. 


2. ADSORPTION MODEL FOR A HOMOGENEOUS SURFACE OF OXIDE 


The potential determining ions H*, and the cations C* and anions A* of the 
basic electrolyte form the following surface complexes: SOH’, SOH;", 
SO'C* and SOH,*A’ where S is the surface metal atom. The concentrations of 
these complexes on the surface are denoted by [SOH’], [SOH,"], 
[SOC] and [SOH,*A’], respectively. [SO] is the surface concentration of the 
free sites (unoccupied surface oxygens). Except for protons, which are located 
in the potential layer yo, the anions and the cations are situated within the layer 
of potential wg. Then we consider the coadsorption of another cation (a strongly 
adsorbed bivalent ion), which is introduced to the system in the form of a MA3 
salt. Its addition causes the formation of complexes [SO'M”’] in the separate M- 
plane (Figure 1). 


Table 1. Surface reactions and equilibrium equations in the FLM 
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According to the mass action law, the equations for equilibrium, and the 
commonly used in literature constants for the reactions occurring at the ox- 
ide/electrolyte interface take the form presented in Table 1, where a, a4, ay and 
ay are the activities of cations, anions, bivalent cations and protons respec- 
tively. 

According to the diagram in Figure 1, the surface charge density, &, must be 
proportional to the sum of the concentrations of the following surface com- 
plexes: 


ô, = B([SOH} ]+[SOH}A~]-[SO~ ]-[SO~C*]-[SO-M”* }) (1) 


The charge of the adsorbed bivalent ions in M-plane has the- following 
meaning: 


5, =2B[SO-M**] (2) 


and the charges of the specifically adsorbed ions of the basic electrolyte in their 
plane, are given by: 


6, = B([SO C*]—[SOH;A ]) (3) 


where B is the surface density in the appropriate units of charge. 
For the whole compact layer, there must be fulfilled the electroneutrality 
condition: 


6, +u +5, +6, =0 (4) 
The value of the diffuse layer charge, ô, is given by: 
5, = B({SO"]—-[SOH; ]-[SO"M** }) (5) 


The total number of the sites capable of forming the surface complexes on the 
surface, N,, is equal to: 


N, =[SO~]+[SOH’]+[SOH? ]+[SO-C*]+[SOH}A™]+[SO"M**] (6) 


The relationships between the potentials and the charges within the individ- 
ual electric layers are following: 


(7) 


5, eas ô, 
IEA N = Wu =AuWo YW, -Wi= 
Ci C3 
where Of is the model parameter. We assume that the potential Wy changes 
linearly with the potential wo. 
Next, we introduce the surface coverages @;'s by the individual surface com- 
plexes (i=0,+,C,A,M) and free sites (i= -), 
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6, =[SOH°J/N, 6, =[SOH>/N, 0, =[SOH;A VN, 
=[SO-C*)/N, 6, =[SO"M* VN, 6_=[SO"J]/N, =1-)'6, - (8) 


The equivalent description of the reactions leading to the adsorption of ions 
onto the free sites SO’ is now following: 


SO” +H* <—>SOH" (9a) 
SO” +2H* <—>SOH; (9b) 
SO” +C* <—>SO°C* (9c) 
SO” +2H* + A~ <—>SOH}A™ (9d) 
SO” +M” <—>SO°-M* (9e) 
And now we have the following equilibrium constants: 
K,; = | ibe Ko = Ke, = . (10) 


int +` print pint Cig int A` prim * prim 
Ki Ka K Ka Kiz K; 


along with the following equilibrium equations: 


KS exp| Se = Ca. 6. (11a) 
kT 8, 

inte A 2ey, | _ (au) 0_ ilb 
ka roofer) -ae i 
sar te} ee (4, )0_ (11c) 
1 rah ey {erated sorsia (y)"(4,)0- at (11d) 

int int 2e Wu (ay )0_ 1 1 
KMKG exp| zeta 1 \- mya (11e) 


The set of the nonlinear eqs 11 can be transformed into another one, having 
the form of multicomponent Langmuir-like adsorption isotherms of ions 8's 
(i=0,+,A,C,M): 
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6, = wh i=0,4,A,C (12) 
DKA 


where f's, (i=0,+,A,C,M), are the following functions of activity of protons, 
bivalent and salt ions in the equilibrium bulk electrolyte, 


fy = exp- re = 23H}, f= fè (13ab) 
ô 2ea 
fe =a spl- an + A, fu =2y exp- Iwl (13cd) 
1 
e eô 
fi =a, erpl- Ye -êe spn) (13e) 


The surface potential function y (pH) defined in the equations for the equi- 
librium of the surface complexation reactions can theoretically be calculated 
from the equation, developed by Bousse and co-workers [15] and used in our 
papers [2-5]. The simplified linear form of this equation is given by, 


_ B 23kT 
Yo = B+1 


which readily shows the difference between the yọ (pH) dependence and the 
Nerstian one. 

Only for relatively large values B, B/(6+1)~1, the potential change corre- 
sponding to one pH unit in this equation becomes that one predicted by the 
Nernst equation i.e. 2.3kT/e Volts for a pH unit. 

The nonlinear equation system 12 can be transformed into a one nonlinear 
equation by taking into account eq 1, 


——(PZC- pH) (14) 


6, = B, Refs thats ES Kulu) jc04,A0.M (15) 
+ Kf, 


This nonlinear equation for & can be solved by means the iteration method, 
to give the value of & for each pH. Having calculated these values, one can 
evaluate the individual adsorption isotherms from eqs 12. 

The Rudzifski-Charmas criterion [3,4] for the common intersection point 
(CIP) can be applied to study the relations between the intrinsic equilibrium 
constants of the reactions in Table 1, and the point of zero charge (PZC) for this 
model. 
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Introducing the notation: PZC = pH, -oy =—logH where H is the activity 
of protons in the bulk solution at PZC, pK'" =-logK™ (i=1,2), 
pK" =-log’K'"™ (i=C,A) and assuming that the concentration of M” is 
very small, the relations reducing the number of the independent equilibrium 
constants take the following form [3,4]: 


PEC =—(pK'" + pK) (16a) 


PEC =—(p'Ke +p'Ki") (16b) 


3. ADSORPTION ON HETEROGENEOUS SURFACES 


It is now commonly realized that the crystallography and chemical compo- 
sition of the actual oxide surfaces do not represent an extrapolation of appro- 
priate bulk crystal properties. The actual (really existing) solid surfaces are 
characterized by a more or less decreased crystallographic order, leading also to 
variations in the local chemical composition. This, in turn, causes variations in 
adsorptive properties of adsorption sites, across the surface. That phenomenon 
known as the "energetic heterogeneity" of the real solid surfaces is believed 
now to be one of the fundamental, common features of the actual solid surfaces 
[16,17]. 

The quantitative measure of the degree of surface heterogeneity in the model 
of the one-site-occupancy adsorption is the differential distribution of the frac- 
tion of surface sites among the corresponding values of adsorption en- 
ergy €, x(€), 


i y(€)de =1 (17) 


where Ae is the physical domain of €. Very often y(e) is approximated by a 
simple analytical function and Ae is assumed to be either (0,+0°) or (—s9,+0°) 
interval for the purpose of mathematical convenience. Such simplifications do 
not usually introduce a significant error in the latter theoretical calculations, 
except for some extreme physical regimes [16]. The exact function (e€) for a 
real physical surface is expected to have a complicated shape in general. How- 
ever, to a first crude approximation, it may be approximated by a simple 
smooth, gaussian-like function. Thus, we will represent y(e) by the following 
gaussian-like function [18], 
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expe} 
[repke] 


centered at €= £", the spread (variance) of which is described by the heteroge 
neity parameter c. (The variance ø is equal to 2c/ v3 ). 

In the case of monolayer adsorption the use of this function leads to the 
Langmuir-Freundlich isotherm which is probably the most commonly applied 
to correlate the experimental data at not very high surface coverages. Using 
such a simple function x(€) is a necessity in practical calculations, to avoid 
introducing many unknown parameters. As it was discussed by Rudziński et al. 
[2] the intrinsic constants K;'s can be written as follows: 


X(E)= (18) 


K,=K. afi, i=0,+,A,C,M (19) 
kT 


where e, is the adsorption (binding) energy of the ith surface complex, and K; 
is related to its molecular partition function. 

The experimentally measured adsorption isotherms have to be related to the 
following averages, 0, , 


6, ({a},T) = f iez f 6,({e},{a},T)x({e}) de,de, de ,dEcdE y (20) 


where {a} is the set of the bulk concentrations {ay, ac, aa}, {e} is the set of the 
adsorption energies {&,€,,€,,&,€y}, Ae is the physical domain of {e },and 
x({e}) is the multidimentional differential distribution, normalized to unity, of 
the number of adsorption sites among various sets { € }. 

In the case of a heterogeneous oxide surface ¢, has different values on dif- 
ferent SO sites of the oxide surface. This is ‘accompanied by the changes in 
K, , constant, but it is generally believed, that these changes are of a secondary 
importance compared to the changes in the chemical bond energy £; . 

Then, for the reasons explained in the previous paper by Rudzinski et al. [4] 
we accept the random model of surface topography, where centres correspond- 
ing to various sets { £ } are randomly scattered on a solid surface. Such random 
centre distribution causes that the probability of finding any other centre close 
to an adsorption centre is the same. As a result a microscopic composition of 
the adsorbed phase close to any centre is the same and it is identical with a 
mean composition of the adsorbed phase on the whole surface { 6,, }. It means 
that all interaction potentials are a function of the averaged concentrations of 
the surface complexes { 8, }. In eq 12 the functions 6,'s have the same form as 


A four layer model for bivalent cations adsorption... 195 


for the homogeneous surface model except for the fact that sets { 8, } should be 
replaced now by sets { 8, }. The adsorption system constitutes a thermodynamic 
entity characterized by one electrostatic capacity. 


El el 


—_—> ĉi 

Figure 2. (A) Schematic visualization of the adsorption energy distributions y,{e,}, 
(i=0,+,A,C,M) for the model assuming high correlations to exist between the adsorption 
energies of various surface complexes. (B) Schematic visualization of the adsorption 
energy distributions y,{¢,}, (i=0,+,A,C,M) for the model assuming lack of correlations 
between the adsorption energies of various surface complexes 


Now, we have to consider another factor, discussed in detail in [5], charac- 
terizing the properties of a heterogeneous surface. This is the degree of correla- 
tion between adsorption energies of various complexes on various centres: 
(1)The case of high correlations between adsorption energies (Figure 2A). 

Though the energies ¢, and ¢,,, change on passing from one centre to an- 

other, their difference A, =€, —&, remains unchanged. Therefore the func- 

tion y({e}) in eq 20 reduces to one-dimensional differential distribution 
xe}. 

(2)The case of lack of correlations between adsorption energies on various 
centres (Figure 2B). 


Then the function y({e}) in eq 20 becomes a product of one-dimensional dis- 
tribution functions. 

The model of a heterogeneous surface taking into account these two extreme 
cases of random surface topography was used in the calculations presented in 
our previous paper [5]. We used it to explain the behavior of adsorption iso- 
therms of bivalent ions. The concentrations of these ions are low, and their 
adsorption characteristics are strongly affected by surface energetic heteroge- 
neity. The logarithm of the adsorbed amount, plotted vs the logarithm of the ion 
concentration in solution is always a Freundlich linear plot with a tangent much 
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smaller than unity. However, at a sufficiently low ion concentration, a transi- 
tion occurs into a Henry's plot, typical for a homogeneous solid surface, with 
the tangent equal to unity. Firstly, we simplified the mathematical calculations 
assuming that Ae is the interval (—s,+ec). Assuming such integration limits 
leads to the isotherm equations which, however, do not reduce to Henry's equa- 
tion for very low surface coverages for model with or without correlations be- 
tween the adsorption energies of various surface complexes. 

However, the equation obtained for the heterogeneity model assuming lack 
of correlations between the adsorption energies - i.e. when are no correlations 
exist between the one-dimensional adsorption energy distributions x,{¢,} given 
by eq 18 - worked well for higher surface coverages giving the Freundlich lin- 
ear plot with a tangent smaller than unity. For such a case we obtained [5] the 
following adsorption isotherm equations for different surface complexes: 


Ya 
it “ate i=0,4+,A,C,M (21) 
where 
Brak; cx aed | (22) 


and where c;'s are the heterogeneity parameters for different surface complexes, 
i=0,+,A,C,M. The way of solving of the equation set 21 was shown in our pre- 
vious paper [5]. 

The general form of adsorption isotherm in the FLM described here is iden- 
tical to eq 21 obtained originally for TLM. The difference appears only in the 
form of the function fy. 

Now, we are going to establish the correlations between the intrinsic equi- 
librium constants in a similar way as in the case of the homogeneous surface 
model, considered in the previous section. The Rudzinski-Charmas criterion 
leads now to the following interrelations [4]: 


kT kT 
2 c, * yin Ce 
se yt os A p n e a it EA (23a) 
oil = kT Ic, $e 
ar Ar 
2 ca yin C 
ET lt Balt g BE Be AT (23b) 
Cy \ Ka Ky c.\ Ka 


In the studies on monolayer adsorption of gases Rudziński and Everett [16] 
showed that a correction for the physically reasonable limited domain Ae, i.e 
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existence of minimum and maximum values £; and £” should be made to ar- 
rive at an isotherm equation reducing correctly to Henry's law. 

Assuming such finite integration limits for the two different models of ran- 
dom surface we obtained [5] two different sets of equations for @,,, corre- 
sponding to two different correlation models. These sets were next used to cal- 
culate the surface coverage by bivalent metal ions at a small ion concentration. 
The surface coverages for other ions have higher values, therefore the equations 
can be used obtained by assuming infinite integration limits. After numerous 
calculations it proved unattainable. The model assuming high correlations be- 
tween the energies of ion adsorption did not reflect existence of two different 
linear parts in the log-log representation and predicted existence of Henry's 
region up to high coverages [5]. The results of our calculations based on the 
model assuming lack of correlation between the energies of various complexes 
predicted a transition from Henry's region (the tangent of log-log plot is equal 
to one) to Freundlich's region (the tangent smaller than one). The problem is 
that the tangent in this transition region reaches values larger than unity [5]. 
Such behavior has never been observed in the experiment. However, the differ- 
ent behavior of the log-log plots suggest, that the degree of the correlations 
between the adsorption energies of ions may affect strongly the behavior of 
these adsorption systems at low ion concentrations. 

Therefore we have decided that the equation describing a real adsorption 
isotherm @,,of the heavy metal ions will be represented by an analytical for- 
mula corresponding to the assumption that x ,(¢;) is the following rectangular 
energy distribution: 


=i,  fore;e<e},e7 > 
X;(E)= ö i (24) 


elsewhere 


The rectangular distribution is a good approximation in the case of strongly 
heterogeneous surfaces [16]. And this is just the case of the bivalent metal ion 
adsorption (Freundlich's plots with kT/c;<<0.9). As the transition from Henry's 
to Freundlich's plot has been observed in adsorption on so different materials, 
we arrive at the following important conclusion. The model assuming high 
correlations between the adsorption energies of various surface complexes is to 
be abandoned in the studies of ion adsorption within the electrical double layer 
formed at the water/oxide interfaces. 

On the contrary, our model calculations corresponding to the assumption 
that no correlations exist between the adsorption energies of different surface 
complexes, could reproduce very well the transition from Henry's to Freund- 
lich's plot [5]. The adsorption isotherm equation obtained by using the rectan- 
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gular distribution function for the random model without correlation takes the 
form [5], 


1+(K°f,)” 1+expy it 
o, s E mit yin r wey =J=M (25) 
1+ (KIS) EFE; L+expyet 


where € has the following form, 


€,.=-kTInK°f,+e°, j=M (26) 


The form of the eq 25 obtained for the TLM [5] will be identical for FLM as 
well. The difference is only in the form of the function fy. 
The others @;, (=0,+,A,C) can be calculated form eq 15 as their coverages 


reach higher values. 
4. RESULTS AND DISCUSSION 


For the FLM model and the heterogeneous surface model without any cor- 
relations between adsorption energies (eq 25), we are going to demonstrate its 
utility to fit the experimental data. Similarly as in our previous investiga- 
tions [5], we took the set of parameters obtained by Davis and Leckie [19] 
analyzing their titrations curves (Table 2). 


Table 7. Values of the parameters obtained by Davis and Leckie [19], who analyzed 
their titration data for the amorphous iron hydroxide in terms of the homogeneous sur- 
face model 


pK‘ =10.70 pK“ =6.90 
PZC=7.90 N,=10 sites/nm? 


c, =1.40 F/m? 1=0.1 mol/dm’° (NaNO;) 
T=25°C ee ee BE co 


We calculated, however, the corresponding pK™ and p’K‘\" values from 
relations 23ab and then adjusted the other parameters: pK, ,¢),¢7', and a, to 
fit the experimental isotherms of Cd”* ions for the three pH values. The results 


of these numerical exercises are shown in Figure 3. 
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Figure 3. Effects of pH on the isotherms of cadmium ions, adsorbed on the amorphous 
iron oxyhydroxide investigated by Benjamin and Leckie [20]. Experimental data for 
three pH values: («) pH=6.6; (@) pH=6.9; (™) pH=7.2. Our theoretical lines (solid) 
were calculated using eq 25 and the same set of the best fit parameters for each pH 


Our theoretical lines were calculated from eq 25 by using for the investi- 
gated system the constants: N, and PZC and the parameters pK", pK" and 
cı collected in Table 2. As to the values of these three parameters their precise 
determination is necessary to analyze the experimental data concerning the 
adsorption of basic electrolyte ions (radiometric curves, potentiometric titration 
curves for different salt concentrations and electrokinetic curves) [4]. And be- 
cause we did not have such data for the studied system, the values of Davis and 
Leckie parameters were used as acceptable approximation because they do not 
affect significantly the log-log adsorption curve of bivalent metal ion adsorp- 
tion. 

Morover, the parameters characterizing the heterogeneous surface model 
without correlations between the energies of adsorption of various surface 
complexes kT/c; (i=0,+,A,C) were assumed to be equal to 0.7, as these values 
do not affect significantly the adsorption of bivalent ions. The values of the 
parameters pK". =10.7 and p’K'"=6.8 were calculated from relations 23ab. 

So, in our numerical exercises we adjused only the three other parameters: 
pK\" =4.4, ei, = -8 kJ/mol, e» =8 kJ/mol. The strategy of application of eq 25 
has been presented in our previous paper [5]. As the slope of the Freudlich part 
of all these isotherms is equal to 0.66, we took the kT/cy as equal to this value 
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to calculate the isotherm in Henry's, transition, and in Freundlich regions for all 
the investigated pH values. 


0.0 0.2 04 0.6 0.8 1.0 1.2 14 


Figure 4. Relations between the correlated parameters pxi* and q,, in the four layer 
model studied. Calculations were made for the system presented in Figure 3, separately 
for the three pH values: (4 ) pH=6.6; (@) pH=6.9; (M) pH=7.2 


To sum up the above considerations, it can be stated, that the log-log ad- 
sorption curves are sensitive to the last three parameters as well as to the pa- 
rameter a,, characteristic for the FLM. a,,, which was accepted as equal to 
one, because our model investigations showed that it is strongly correlated with 
the adsorption constant pKi". The extent of this correlation is given in Fig- 
ure 4. 
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